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Synchronized Magnetization Oscillations in F/N/F Nanopillars 

Kiwamu Kudo*, Rie Sato and Koichi Mizushima 

Corporate Research and Development Center, Toshiba Corporation, Kawasaki, 212-8582, Japan 

Current-induced magnetization dynamics in a trilayer structure composed of two ferromagnetic free layers and a 
nonmagnetic spacer is examined. Both free layers are treated as a monodomain magnetic body with an uniform 
magnetization. The dynamics of the two magnetizations is modeled by modified Landau-Lifshitz-Gilbert equations 
with spin-transfer torque terms. By solving the equations simultaneously, we discuss their various solutions in 
detail. We show that there exists the synchronous motion of two magnetizations among the various solutions; 
the magnetizations are resonantly coupled via spin-transfer torques and perform precessional motions with the 
same period. The condition to excite the synchronous motion depends on the difference between the intrinsic 
frequencies of the two ferromagnetic free layers as well as the magnitude of current. 

KEYWORDS: spin-transfer torque, current-induced magnetization dynamics, Landau-Lifshitz-Gilbert equation, synchro- 
nization, magnetization oscillations 



1. Introduction 

Current-induced magnetization excitations in mag- 
netic nanopillars, predicted by Slonczewski 1 and Berger 2 
in 1996, have been extensively studied. The basic 
structure of nanopillars is ferromagnetic/nonmagnetic/ 
ferromagnetic (F/N/F) trilayer structure. It has been 
demonstrated experimentally that a spin-polarized dc 
current through a "free" ferromagnetic layer can reverse 
its magnetization. 3-8 Moreover, recent experiments have 
shown that a coherent precession of the magnetization 
at GHz frequencies is induced by a spin-polarized cur- 
rent and that its frequency of precession depends on 
applied field and current density. 9-11 These magnetiza- 
tion dynamics, current-induced magnetization switching 
(CIMS) and coherent precession, provides the possibility 
to utilize the spin-transfer phenomena for various appli- 
cations, such as magnetic random access memory cells, 
nanometer-sized microwave generators, and so on. 

Most experimental studies about CIMS and microwave 
generation have concerned an "asymmetric" trilayer con- 
sisting of a free thin ferromagnetic layer and a thick ferro- 
magnetic layer with a fixed magnetization. Accordingly, 
many theoretical works also have concerned the asym- 
metric structure. 12-17 By active experimental and the- 
oretical studies for the past ten years, many properties 
of magnetization dynamics in the asymmetric structure 
have already been clarified. On the other hand, the mag- 
netization dynamics in a "symmetric" structure where 
both ferromagnetic layers play the same role as free lay- 
ers has received little attention. It is reported recently 
that precessional dynamics of the thick layer as well as 
the thin layer are excited by current in the usual asym- 
metric sample. 18 Therefore, it is important to examine 
what kind of magnetization dynamics can be excited in 
the two-free-layers structures. 

In this paper, we examine current-induced magnetiza- 
tion dynamics in the F1/N/F2 nanopillars where both 
Fl and F2 are ferromagnetic free layers and N is a non- 
magnetic spacer as illustrated in Fig. 1. By the current 
passing through the trilayer, the magnetizations of Fl 
and F2, M\ and M2, interact with each other via spin- 
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Fig. 1. Schematic of the F1/N/F2 trilayer. Fl and F2 are thin 
ferromagnetic "free" layers. N is a nonmagnetic spacer. J is cur- 
rent density perpendicular to plane. Fl and F2 are treated as a 
monodomain fcrromagnet. 



transfer torques, 1 and perform various motions. We dis- 
cuss these motions in detail comparing with the magne- 
tization dynamics in an asymmetric structure. We show 
that among these various motions there exists a syn- 
chronous motion of two magnetizations. In the motion, 
the two magnetizations perform a stable precessional mo- 
tion with the same period, analogous to the synchronized 
oscillation of two coupled nonlinear oscillators. The syn- 
chronization of magnetizations may be efficient for rais- 
ing the power of microwave radiation emitted by the 
magnetic multilayers which function as microwave gen- 
erators. 

2. Model 

With regard to two ferromagnetic layers in the Fl/ 
N/F2 trilayer, we assume that both Fl and F2 have an 
uniaxial anisotropy in the z direction. We assume further 
that Fl and F2 have almost identical properties except 
that there is a difference between the magnitude of the 
effective magnetic fields acting on Mi and A^- Within 
monodomain approximation, we describe the dynamics 
of the two magnetization vectors, M\ and M2, by the 
modified Landau-Lifshitz-Gilbert equations 



dm, 



drrij 



(1) 



— — = -m, x hQ + am l x — — 

d,T clt 

with the spin-transfer torque term I\ . Here we have used 
the dimensionless form of LLG equations. In eqs. (1), 
rrii = Mi/ M s , M s is the saturation magnetization, times 
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r is measured in units of ( i yH^)~ 1 , is the magni- 
tude of the uniaxial anisotropy field acting on Mi, 7 
is the gyromagnetic ratio, and a is the dimensionlcss 
Gilbert damping constant of which we always take 0.01 in 
this paper, h^g is the effective magnetic field normalized 
by Hu and given by 

h® = + h^m iz )z - hfm lx x, (2) 

where ftg Xt is the external field applied along the z 
direction, b$ is the in-plane uniaxial-anisotropy field 
{h^ = 1.0), and h£\~ AttMs/H^) is the effective out- 
of-plane anisotropy field due to the film shape of ferro- 
magnets. 

The spin-transfer torque term Ti acting on Mi can be 
written in the form 1, 13 

r, = —ajrrii x (mi x m 2 ). (3) 

The parameter aj represents the strength of the spin- 
transfer torque and is proportional to current density 
J. We assign the positive value of aj for the current 
flowing from Fl to F2. It is noticed that the strength 
of the spin-transfer torque acting on mi and m 2 are 
equal because we have assumed that Fl and F2 have 
almost identical properties. By means of eqs. (3), the two 
magnetizations can interact with each other. The motion 
of m 2 contributes to that of mi via Ti, and vice versa. 
Due to the form of eqs. (3), I\ = when mi and m 2 
are parallel. 

We solve eqs. (1) numerically and examine the current- 
induced dynamics of m^ on the unit sphere \mi\ = l. 19 
We conduct the calculation using mainly the two param- 
eter sets for the effective fields; (i) ftg Xt = ^cxt = 2.0, 
h ( u 1] = 1.0, h ( u } = 1.0 + Aft u , h { p ] = ft p 2) = 0.0, and 
(ii) fc£> = 1.0, = 1.0 + Aft ext , h ( u 1} = h u 2) = 1.0, 
hp^ = h[P = 10.0. In the parameter set (i), we ne- 
glect the effect of out-of-plane demagnetizing fields ftp ' 
for simplicity. On the other hand, the parameter set (ii) 
is realistic and experimentally realizable. In the both 
sets of parameters, we model the difference between the 
magnitude of effective magnetic fields acting on the two 
magnetizations by Aft u and Aft oxt , respectively. We as- 
sume that Ah u > and Aft oxt > 0. Since Ah u > or 
A/iext > 0, we consider the trilayers where the magnitude 
of effective magnetic field acting on m 2 is larger than 
that acting on mi in both models (i) and (ii). As we dis- 
cuss below, the difference of the effective magnetic fields, 
such as Ah u and Ah ex t, has the essential role to arise a 
synchronous motion of the two magnetizations. By using 
the simpler model (i) , we clarify the characteristics of the 
synchronous motion of the magnetizations. By using the 
model (ii), we show that the synchronous motion exists 
even in a realistic setup. The parameters in (ii) cover the 
essential features of magnetization dynamics in the al- 
most symmetric F/N/F trilayers including the two free 

(i) (i) (i) (i) 

layers with h u < h e ^ t <C hp . Note that ftp is several 
tens of times larger than ftu ' in usual film ferromagnets. 
Moreover, Aft cxt can be realized by a exchange-bias field 
as that in spin-valve pillars. 



In both models (i) and (ii), the equilibrium direction of 

(i) 

the two magnetizations along h c g corresponds to m^ = z 
in the absence of current. Accordingly, it would be ap- 
propriate that mi — z as the initial value at r = 0. 
However, when mi and m 2 arc completely parallel at 
r = 0, the spin-transfer torque, eq. (3), does not work 
at all and any magnetization dynamics are not excited. 
Therefore, we permit that the two magnetizations shift 
slightly from their equilibrium directions and arc initially 
not parallel each other. This is justified by taking into 
account the effect of finite temperature. In the real ma- 
terials with temperature T 7^ 0, thermal fluctuations 
always exist and then magnetizations fluctuate around 
their equilibrium points all the time. Therefore, we as- 
sume that two magnetizations are initially in the vicinity 
of their equilibrium points and are not parallel, although 
we do not take account of thermal fluctuations when we 
solve eqs. (1). 

3. Synchronized Precession of Two Magnetiza- 
tions 

In this section, we discuss the synchronous motion of 
the magnetizations by showing the results of the calcu- 
lation in model (i), i.e., a = 0.01, ft^j. = ft^ = 2.0, 
hu 1] = 1.0, h u 2) = 1.0 + Aft u , ft p x) = ft p 2) = 0.0. We con- 
sider only the case that the applied current is positive, 
aj > 0, because the synchronization of two magnetiza- 
tions do not occur when aj < in model (i). 

3.1 Phase diagram 

Depending on aj and Aft u , several distinct types of dy- 
namical modes as the solutions of eqs. (1) are excited by 
the spin-transfer torque. Figure 2 shows the dependence 
of steady state solutions on aj and Aft u . In the region 
labeled S, the synchronized precession of two magnetiza- 
tions on which we will focus in this section occurs. In the 
other regions which are labeled by P and W, any kind of 
coherent precessions do not exists. The region P denotes 
the region where a static parallel configuration of two 
magnetization occurs. The region W denotes the region 
where two magnetizations behave chaotic. 

Figure 3 shows the typical behaviors of mi and m 2 
in the region W. It is found that the behavior like mag- 
netization switching arises irregularly. To obtain Fig. 3, 
we have used m lz (0) = 0.998, m 2z (0) = 0.999, and 
(p(0) = -tt/6 as initial conditions. Here, ip = y 2 — <Pi 
where ipi is the azimuthal angle of m^. In other words, ip 
is the angle between m^ and m 2 which are the projec- 
tions of mi and m 2 onto the xy-plane, respectively; see 
Fig. 5(d). In the horizontal axis of Fig. 3, /o is the intrin- 
sic frequency of mi in the absence of current, /o = loq/ 
2tt, where uj is given by the well-known formula, 

^0 = v / (^ x ) t + ^ 1) )(^ x ) t + ^ 1) + 4 1) )- (4) 

u>o can be determined by the ferromagnetic resonance 
(FMR) experiment. We have introduced /o as the stan- 
dard precessional frequency. In model (i), /o = 3.0/27T. 

In the region S, both mi and m 2 precess around the 
z axis with circular forms. The circular trajectories of 
the precessional motion of mi and m 2 are shown in 
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Fig. 2. (Color) Dynamical stability diagram for aj and Ah u . S 
denotes synchronized precession of two magnetization, P paral- 
lel configuration, and W chaotic motions. The numerical values 
in the legend represent the steady-state values of m\ z . X repre- 
sents parallel configuration or chaotic motions. The broken line 
L represents the threshold for current-driven excitations. 
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Fig. 3. Typical behaviors of mi and m.2 in the region W. Time 
evolution of mi z (j) is shown. The parameters are chosen as 
aj = 0.08 and A/i u = 0.8. The initial conditions are chosen 
as m lz (0) = 0.998, m 2z (0) = 0.999, and ip = tt/6. 



Fig. 4(a). The steady precessions occurs after transient 
behaviors. The behaviors are shown in Fig. 4(b) where 
the time evolutions of m\ z (r) and mi z (j) are plotted. 
Figure. 4(c) shows the time evolutions of the x compo- 
nents, mia;(r) and TO2 X (t)) in the steady state at inter- 
vals between for = 300 and 305. As can be recognized 
from Fig. 4(c), the precessional frequencies of mi and 
mi are equal. In other words, the precessional phase is 
mutually locked. That is why we name the magnetiza- 
tion dynamics in the region S "synchronized precession" 
of two magnetizations. 

3. 2 Synchronized precessions 

We discuss the properties of the synchronized preces- 
sions in detail. 

The precessional amplitudes become large as aj and 
Ah u are increased. The dependence of the amplitude of 
mi on a j and Ah u has been shown in Fig. 2. The legend 
in Fig. 2 represents the values of m\ z in the synchronized 
precessional state. Although it is not illustrated, the am- 
plitude of m.2 has the similar dependence on aj and Ah u 
as that of mi. The difference between the precessional 
amplitude of mi and m2 is that the amplitude of m2 is 
larger than that of mi as is shown in Fig. 4(a) for any 
a j and Ah u . This property that m\ z > m,2 Z is essential 
for the stability of the synchronous precession. 

In the synchronized precession, the two magnetizations 
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Fig. 4. Typical behaviors of mi and m.2 in the region S. The 
parameters are chosen as aj = 0.06 and A/i u = 0.4. The initial 
conditions are chosen as mi z (0) = 0.998, m2z(0) = 0.999, and 
ip = 7r/6. (a) Bird's-eye view of trajectories in the steady state, 
(b) Time evolution of rm z (r). (c) Time evolution of m^fr) at 
intervals between for = 300 and 305. 



precess around the z axis with the same period. There- 
fore, the phase shift ip is locked. The magnitude of the 
locking phase shift depends on a j and Ah u as is shown in 
Fig. 5(a). The numerical values in the legend of Fig. 5(a) 
represent the magnitude ofip/n. The phase shift becomes 
small as the magnitude of current increases. 

The angle between mi and m-2, <fi — arccos(mi • m^), 
also depends on aj and Ah u in the synchronized preces- 
sional state. The dependence is shown in Fig. 5(b). The 
numerical values in the legend of Fig. 5(b) represent the 
magnitude of 4>/w. The angle between the two magne- 
tization has the largest value in the precessional region 
where low aj and high Ah u . 

The precessional frequency w s in the synchronized pre- 
cession is tunable by current. In Fig. 5(c), we have shown 
the dependence of tu s on aj and Ah u . The legend in 
the figure represents the value of uj s . It is found that 
u s monotonously reduces as the magnitude of current 
increases. 

Next, let us discuss the conditions for the excitation 
of the synchronized precession. 

In Fig. 2, we have shown the threshold for the current- 
driven synchronized precession. The threshold is denoted 
by the broken line L. We have found that L can be esti- 
mated by aj — ah]: where h^ is given by 

%=hW+hg\ (5) 

Therefore, it is necessary that a.j > ah^(— a(3.0+ Ah u j) 
for the synchronized precession. It is noticed that the 
condition coincides with that for magnetization excita- 
tions in an asymmetric structure which can be obtained 
by a modified LLG equation; see eq. (17) in ref. 12. L 
corresponds to the threshold for magnetization excita- 
tions of m-2 in a m-i-pinned asymmetric trilayers. The 
reason that the current threshold L depends on the ef- 
fective magnetic fields acting on tti2 can be understood 
as follows; when the two magnetizations, mi and m-2, are 
initially almost parallel and a positive current (aj > 0) is 
applied, the direction of the spin-transfer torque is such 
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Fig. 5. (Color) (a) Dependence of the phase shift ip on aj and 
Ah u - (b) Dependence of the angle between the two magneti- 
zations <f> on a j and Ah u . (c) Dependence of the precessional 
frequency lu b in the synchronized precession on aj and Ah u . (d) 
Definition of ip and (p. 

that it induces the dynamics of m-i at first. On the other 
hand, when a negative current (aj < 0) is applied, the di- 
rection of the spin-transfer torque is such that it induces 
the dynamics of mi at first. 

It is noticed that the region P has spread in all ranges 
of aj for very small Ah u in Fig. 2. Therefore, it is nec- 
essary that Ah u ^ as well as that aj > ah]: for the 
excitation of the synchronized precession. In other words, 
it is necessary for the synchronized precession that the ef- 
fective magnetic fields acting on the two magnetizations 
are different. We generalize this necessary condition in 
§4.3. 

Finally, we compare the results discussed above with 
the ones obtained by the macrospin model in the asym- 
metric structures. In the modified LLG model for the 
asymmetric structures, coherent magnetization preces- 
sions are possible only when h u < h ey ± < h u + h p . 17 
Here /i e xt, h u (= 1), and h p are the dimensionless exter- 
nal magnetic field parallel to planes, in-plane uniaxial 
anisotropy field, and out-of-plane anisotropy field which 
act on the magnetization of the free layer, respectively. 
Accordingly, the out-of-plane anisotropy field h p is im- 
portant for current-induced magnetization precessions in 
the asymmetric structure. On the other hand, the syn- 
chronized magnetization precessions in the symmetric 
structure are possible even in the absence of hp . Ac- 
tually, we have neglected the effect of h p in this section. 

4. Magnetization Dynamics in Realistic Trilay- 
ers 

In this section, we discuss the magnetization dynam- 
ics including synchronized precessions in realistic trilayer 
structures. As a realistic setup, we use the set of param- 
eters (ii): a = 0.01, = 1.0, = 1.0 + Ah cxt , 
h ( u } = = 1.0, and h p 1} = h p 2) = 10.0. 



CP LI L2 




Fig. 6. (Color) Dynamical stability diagram for aj and Ah CX f 
S denotes synchronized precessions, Q quasi-periodic motions, 
OPP out-of-precessions, CP large-angle clamshell precessions, W 
chaotic behaviors, and P parallel configuration. 20 The numerical 
values in the legend represent the values of (mi z ). X represents 
parallel configuration or chaotic behaviors. The broken lines, LI 
and L2, represent the threshold for current-driven excitations. 

4-1 Phase diagram 

In Fig. 6, the dynamical stability diagram for aj 
and A/i cxt is shown. S, Q, OPP, CP, P, and W repre- 
sent the region where synchronized precessions, quasi- 
periodic motions, out-of-plane precessions, clamshell pre- 
cessions, parallel configuration, and chaotic behaviors oc- 
cur, respectively. The numerical values in the legend rep- 
resent the values of (mi z ) in the synchronized precession, 
where (■ • • ) is the one cycle average. The broken lines, LI 
and L2, are the threshold for current-driven excitations. 
LI and L2 can be estimated by aj = —ah]: 1 where 

/* C L1 = a£ + + = 7 -°< (e) 

and a j = ah]: 2 where 

= + h{? + *L = 7.0 + A/w, (7) 

respectively. It is then necessary that aj < —ah]? 1 or 
aj > ah]: 2 for magnetization excitations. 

Let us first discuss the positive current region (aj > 0) 
in Fig. 6. 

In the region S, both mi and reach the steady 
precessional states after transient behaviors as shown in 
Fig. 7(a) and 7(b). In the steady state, the two mag- 
netizations perform precessions around the z axis with 
elliptic forms, whose trajectories are shown in Fig. 7(c). 
The elliptic forms of the precessions result from the ef- 
fect of the effective demagnetizing field hp . As is shown 
in Fig. 7(d), the precessional period of mi and m.2 are 
identical. That is, the phase shift between mi and 
is locked. Therefore, the synchronized precessions exist 
even in a realistic setup. The precessional amplitudes be- 
come large as aj and A/i oxt are increased. This tendency 
is the same as that discussed in the previous section. 

For larger aj or A/i oxt , the phase locking become 
weak and quasi-periodic motions occur. The region where 
the quasi-periodic motions occur is denoted as Q in 
Fig. 6. The typical motions of mi and m,2 in the re- 
gion Q are shown in Fig. 8. Both mi and m^ per- 
form precessions around the z axis with several char- 
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Fig. 9. (Color) (a) Typical power spectrum in the region S. The 
parameters are chosen as (aj, Ah ex t) = (0.24813,0.6). (b) Typ- 
ical power spectrum in the region Q. The parameters are chosen 
as (oj,Aftext) = (0.24814,0.6). 



Fig. 7. Typical behaviors of mj and in the region S. The 
parameters are chosen as aj = 0.2 and A/i cx t = 0.6. The ini- 
tial conditions are chosen as mi z (0) = 0.999, m,2 z (0) = 0.998, 
and ip = 7r/3. (a) Time evolution of m\ z {r). (b) Time evolution 
of Tnizir). (c) Bird's-eye view of the trajectories in the syn- 
chronized precession, (d) Time evolution of mi z (r) at intervals 
between a/or = 1.2 and 1.24. 
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Fig. 8. Bird's-eye view of the trajectories in the region Q. The 
parameters are chosen as aj = 0.3 and Ah ex t = 0.6. 

acteristic frequencies. The typical power spectrum for 
the quasi-periodic motions is shown in Fig. 9(b). I(u>) = 
lim-T^oo ^ (m* x (u>)rni x (u))) is plotted, which is obtained 
by the Fourier transform of mi x {r) in the steady state. 
It is found that several characteristic frequencies coexist. 

The phase boundary between the region S and 
Q is very sharp. In Fig. 9(a) and 9(b), we have 
shown the power spectrum for the steady-state mo- 
tions at (aj,Ah cxt ) = (0.24813,0.6) and (aj,Ah cxt ) = 
(0.24814,0.6), respectively. It is determined by the very 
slight difference in aj which motions will occur be- 
tween synchronized precessions and quasi-periodic mo- 
tions. The change to quasi-periodic motions from syn- 
chronized precessions is not continuous for aj, i.e., the 
change is drastic. 

In the region OPP, m-2 performs the out-of-plane pre- 
cession. On the other hand, m\ oscillates only around the 
vicinity of the positive z direction. Those behaviors are 
shown in Fig. 10(a). As a matter of fact, like the mag- 
netization dynamics in the region S, m\ and mi also 
perform precessions with the same period in the region 
OPP. However, because one of them has a small pre- 
cessional amplitude, the magnetization dynamics of the 
whole system is almost an out-of-plane precession. 

For too large aj, chaotic behaviors of mj and ra-i oc- 
cur. The typical behaviors are shown in Fig. 11(a) and 
11(b) where time evolutions of m\ z {j) and m2 Z (r) are 
plotted, respectively. m\ randomly changes its in-plane 




Fig. 10. Bird's-eye view of the trajectories in the region OPP. 
(a) Out-of-plane precession of rri2- The parameters are chosen 
as aj = 0.2(> 0) and A/i ox t = 1.2. (b) Out-of-plane precession 
of mi. The parameters are chosen as aj = — 0.45(< 0) and 
A/i cxt = 0.6. 
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Fig. 11. Typical behaviors of mi z (r) in the region W. The ini- 
tial conditions are chosen as mi z (0) = 0.999, m2 z (0) = 0.998, 
and tp = n/3. (a)-(b) Typical behaviors in the region W where 
aj > 0. The parameters are chosen as aj = 0.4 and A/i ex t = 0.6. 
(c)-(d) Typical behaviors in the region W where aj < 0. The pa- 
rameters are chosen as aj = —0.35 and A?t ex t = 0.6. 

precessional amplitude. On the other hand, m-i hops be- 
tween out-of-plane and in-plane trajectories irregularly. 

Next, let us discuss the negative current region (aj < 
0) in Fig. 6. 

In the region P, due to the smallness of current, any 
kind of magnetization dynamics is not excited. Just the 
static parallel configuration, mi = m.2 = z, occurs in 
the steady state. 

For the current exceeding the threshold LI, mi per- 
forms an in-plane large-angle clamshell precession and, 
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Fig. 12. Large-angle clamshell precession of mi which occurs in 
the region CP. The parameters are chosen as aj = —0.085 and 
Ah cxt = 0.8. 
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Fig. 13. Dependence of the processional frequency u>i on aj when 
Ah cx t = 0.6. The broken line represents u>o — 4.9 determined by 
eq. (4). 



on the other hand, m 2 moves only in the vicinity of the 
z axis. These typical behaviors are shown in Fig. 12. The 
behaviors become unstable at aj = —0.1. 

For aj < —0.1, out-of-plane precessions of mi or 
chaotic behaviors occur. In the region OPP, mi per- 
forms the out-of-plane precession. On the other hand, 
m.2 oscillates only around the vicinity of the positive z 
direction. Those behaviors are shown in Fig. 10(b). In 
the region W, chaotic behaviors of mi and m 2 occur. 
The typical behaviors are shown in Fig. 11(c) and 11(d) 
where time evolutions of m\ z (r) and m 2z (t) are plotted, 
respectively, m^ hops between out-of-plane and in-plane 
trajectory irregularly. m 2 randomly changes its in-plane 
precessional amplitude. It is noted that according to the 
sign of current, the behavior of mi and m 2 is contrastive 
as can be recognized by comparing Fig. 10(a) with 10(b) 
and comparing Fig. ll(a)-(b) with ll(c)-(d). These be- 
haviors in the region OPP and W reflect that the positive 
current region is the m.2-excited region and the negative 
current region the mi-excited region. 

In the negative current region (a,/ < 0), the magne- 
tization dynamics of the whole system is almost deter- 
mined bymj because m 2 oscillates with very small am- 
plitudes as mentioned above. Moreover, the obtained co- 
herent behaviors in aj < are similar to those obtained 
by a macrospin model in an asymmetric structure: in- 
plane precessions and out-of-plane precession. 9 Accord- 
ingly, roughly speaking, the magnetization dynamics is 
very similar to those in an asymmetric structure, espe- 
cially its in-plane high magnetic field region; see Fig. 3 
in ref. 9, for example. The difference is that there exists 
chaotic behaviors in our model. 

4-2 Microwave frequency 

Figure 13 shows the dependence of lc\ on current aj. 
oji is the precessional frequency of mi. Aft, e xt is fixed as 
A/iext = 0.6. The dependence of io\ on aj reflects the 
behaviors shown in Fig. 6 with A/i cxt = 0.6. 

For the negative current regime (aj < 0), we have ob- 
tained the dependence of microwave frequency on current 
similar to the one which can be obtained in the asym- 
metric structure. 9 In the in-planc precessional region, the 
precessional frequency uj\ decreases as the magnitude of 
current \aj\ increases. In the out-of-plane precessional re- 
gion, uoi increases as \a.j\ increases. wi can not be defined 
in the region W because the characteristic frequency does 



not exist in chaotic behaviors. Therefore, in Fig. 13, we 
have not plotted the data in the region W or have labeled 
several points by cross marks. 

For the positive current regime (aj > 0), we have 
plotted the dependence of ioi on aj in the synchro- 
nized precession. The synchronized precessional fre- 
quency monotonously decreases as the magnitude of cur- 
rent is increased. This tendency is the same as that dis- 
cussed in the previous section; see Fig. 5(c). 

4-3 Conditions for the excitation of synchronized pre- 
cessions 

As can be found from Fig. 6, it is necessary that 
A/i C xt 7^ as well as aj > ah]: 2 for the excitation of syn- 
chronized precessions. That is, the difference between the 
effective fields acting on two magnetizations is needed. 
This condition can be generalized; it is necessary that 
there exists the deviation between the intrinsic frequency 
of Fl and F2. The intrinsic frequency of rrii is given 

by 4° = V / (^ ) t + ^u ) )(^ ) t + ^ ) +4 4) )- Accordingly, 
the above condition can be written 
have checked the condition using several sets of param- 
eters. The condition that lu^ ^ lu^ can be understood 
in the following way. If Wq 1 ^ = to^\ Mi and M 2 tend to 
arrange their motions and become parallel to each other. 
As the result, the spin-transfer torque I\ can not work. 
This condition is always applicable when Mi and M 2 
have their equilibrium points in the same direction in 
the absence of current. 

It is noted that synchronized precessions can not ex- 
ist when the deviation of intrinsic frequencies, Aw = 
cjq 2 ^ — is too large; Aw ^> 0. Aw originates in the 
deviation of the effective magnetic fields such as Ah u or 
A/i ext . As it is recognized in Fig. 2 or 6, Alo must be in a 
limited range, otherwise chaotic behaviors occur. One of 
the reason is that the phase locking of two magnetization 
precessions is impossible when Aw is too large. 

4-4 Comparison with macrospin models in an asymmet- 
ric structure 

We compare the behaviors obtained by eqs. (1) with 
the ones obtained by the macrospin model in an asym- 
metric structure. It is known that two types of coherent 
magnetization motions of a free layer can be obtained 
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by an modified LLG equation in an asymmetric struc- 
ture: in-plane precessions and out-of-plane precessions. 
Chaotic behaviors like telegraph noise can not be ob- 
tained by such an macrospin model because the degree 
of freedom of the magnetization is only two. Accordingly, 
the "W" phase experimentally observed by Kiselev et al. 9 
can not be obtained. In our model, however, chaotic be- 
haviors can be obtained by eqs. (1) since the degree of 
freedom of the system is four. The region W appearing 
the negative current regime in Fig. 6 may corresponds to 
the "W" phase discussed in ref. 9. 

5. Conclusions 

Current-induced magnetization dynamics in a sym- 
metric trilayer structure consisting of two ferromagnetic 
free layers and a nonmagnetic spacer is examined. We 
have treated the two free layers as a monodomain fer- 
romagnet and calculated the magnetization dynamics by 
modified LLG equations. We have found that various be- 
haviors of the two magnetizations arise depending on the 
applied current and the deviation of effective magnetic 
fields acting on them. Especially, there exists the syn- 
chronized precessions of two magnetizations among the 
various behaviors. For the excitation of the synchroniza- 
tion, current must exceed the threshold whose magni- 
tude is the same as that for magnetization excitations 
in an asymmetric structure. It is also necessary that the 
intrinsic frequencies of the two magnetizations are differ- 



ent; Wq 1 ' ^ u)q'' . Moreover, the deviation of frequencies, 
Alu = Wq 2 ' — lo^\ must be in a limited range for the 
phase locking of two magnetization precessions. Utiliz- 
ing the synchronization in a symmetric trilayer may be 
one of the possible ways to raise the microwave emission 
power of the spin-transfer oscillators. 21-23 
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